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Abstract. We discuss the microwave second-harmonic generation in high-density bulk MgB2, prepared by 
the reactive liquid Mg infiltration technology. The intensity of the harmonic signal has been investigated 
as a function of temperature and amplitudes of the DC and microwave magnetic fields. The results are 
discussed in the framework of a phenomenological theory, based on the two-fluid model, which assumes 
that both the microwave and static magnetic fields, penetrating in the surface layers of the sample, weakly 
perturb the partial concentrations of the normal and superconducting fluids. We show that, in order to 
account for the experimental results, it is essential to suppose that in MgB2 the densities of the normal 
and condensed fluids linearly depend on the temperature. 

PACS. 74.20.De Phenomenological theories (two-fluid, Ginzburg-Landau, etc.) - 74.25.Nf Response to 
electromagnetic flelds (nuclear magnetic resonance, surface impedance, etc.) - 74. 70. Ad Metals; alloys and 
binary compounds (including A15, MgB2, etc.) 



1 Introduction 

Investigation of the nonlinear microwave response of su- 
perconductors is of great interest for both fundamental 
and applicative aspects. From the basic point of view, it 
allows determining the electron-phonon scattering time TQj 
12 , the temperature dependence of the condensed and nor- 
mal fluids [^^11^ and checking the presence of weak links in 
the sample ln]|^]|7||g) . From the technological point of view, 
it is well known that the occurrence of nonlinear effects at 
high input-power levels limits the application of supercon- 
ductors in passive microwave (mw) devices, while it has 
important implications in active devices |MTU| . For these 
reasons, it is of importance estimating the mw current 
and field intensities at which nonlinear effects are signifi- 
cant and recognizing the mechanisms responsible for the 
nonlinearity at high frequencies. 

Since the discovery of superconductivity in magnesium 
diboride, the scientific world has devoted large attention 
to this superconductor because of its peculiar properties, 
related to the two-gap structure of the electronic states, as 
well as its potential in technological applications ' 11II12I 
ITS] . It has been shown that the presence of the two gaps af- 
fects several properties of the MgB2 superconductor, such 
as the fluxon structure [T^HT^ and the temperature de- 
pendencies of the specific heat QHIj critical fields [TTItT^ . 
condensed-fluid density T9"'20'. The advantages of using 
MgB2 for technological applications are the simple crys- 
tal structure, the malleability and ductility, due to the 
metallic nature, the relatively high critical temperature. 



the large coherence length, which makes the materials less 
susceptible to structural defects like grain boundaries |lll 
I21II22| . One of the main effects of the nonlinear electro- 
magnetic response of superconductors is the emission of 
signals at harmonic frequencies of the driving field j23 |. 
Previous studies carried out in some bulk MgB2 samples 
have shown that the microwave harmonic emission is sig- 
nificant in the whole range of temperatures below Tc and 
exhibits an enhanced peak at temperatures close to 
24.25 . It has been shown that the harmonic emission at 
low temperatures is related to processes occurring in weak 
links, while a different mechanism is responsible for the 
near-Tc peak |24| . Furthermore, measurements of second- 
harmonic (SH) signal at DC magnetic fields smaller than 
the lower critical field, have shown that the SH signal due 
to weak links strongly depends on the specific properties 
of the sample ^11^ . 

Recently, we have investigated the mw SH emission in 
three bulk samples of MgB2 obtained by reactive infiltra- 
tion of liquid Mg in powdered B preforms |27| , which are 
characterized by different mean size of the grains (100, 40, 
1 /im). We have shown that the intensity of the low-field 
and low-T SH signal strongly depends on the size of grains 
and, in particular, that the sample with the smallest grain 
size does not exhibit detectable low-T SH signal, prov- 
ing that in such sample the nonlinear processes in weak 
links are not enough effective in the mw SH generation 
|26j . In this paper, we investigate the SH emission at tem- 
peratures close to Tc in the sample with grain mean size 
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of « 1 /im. The investigation is carried out at low mag- 
netic fields, where nonlinear effects arising from motion 
of Abrikosov fluxons do not play a significant role in the 
harmonic emission. The results are discussed in the frame- 
work of a phenomenological theory, based on the two-fluid 
model, which assumes that both the mw and static mag- 
netic fields, which penetrate in the surface layers of the 
sample, weakly perturb the partial concentrations of the 
normal and superconducting fluids [310] . We show that the 
results expected from this model satisfactorily account for 
the experimental data if a linear temperature dependence 
of the normal and condensed fluid densities is supposed. 



2 Experimental Apparatus and Sample 



3 Experimental Results 

The SH emission has been investigated as a function of the 
temperature, DC magnetic field and input power level. In 
order to disregard SH signals due to nonlinear fluxon dy- 
namics, the attention has been devoted to the SH response 
at low external fields. Before any measurement was per- 
formed, the sample was zero-field cooled (ZFC) down to 
a desired temperature value and then the magnetic field 
was applied. 

Fig. 1 shows the SH signal intensity as a function of 
the temperature in the MgB2 sample, at different values 
of the DC magnetic field. The SH emission is visible, with 
our experimental sensibility, only in a restrict range of 
temperatures below Tc (the dashed line indicates the noise 
level). The SH-vs.-T curves exhibit a peak just below 
on increasing Hq, the peak broadens and its maximum 
shifts towards lower temperatures. 



The investigated MgB2 sample was prepared by reactive 
liquid Mg infiltration technology ^TP. Micron size amor- 
phous boron powder (Grade I, 98% purity. Stark H.C., 
Germany) and pure liquid magnesium are inserted in a 
stainless steel container, with a weight ratio Mg/B over 
the stoichiometric value (~ 0.55); the container was sealed 
by conventional tungsten inert gas welding procedure, with 
some air trapped inside the B powder; a thermal anneal- 
ing at 900 °C for 30 min was performed. It has been 
shown that MgB2 samples produced in this way have a 
density of 2.40 g/cm'^ and consists predominantly of 1- 
micron size grains, although larger grains, of a few mi- 
crons in size, are also present (28) . From the final product, 
we have extracted a sample of approximate dimensions 
2 X 3 X 0.5 mm"^, whose largest faces have been mechani- 
cally polished to obtain very smooth surfaces. 

The sample has been previously characterized by mea- 
suring the AC susceptibility at 100 kHz. From these mea- 
surements we have found that the transition can be de- 
scribed by a Gaussian distribution function of Tc, cen- 
tered at Tco = 38.5 K with ax^ = 0.2 K. Furthermore, 
we have determined the lower and upper critical fields; 
in the range of temperatures of about 5 K below Tc, the 
critical fields show linear temperature dependence with 
dHc2/dT « 1.5 kOe/K and dHd/dT w 6.5 Oe/K. 

In order to detect the SH signal, the sample is placed in 
a bimodal cavity, resonating at the two angular frequencies 
to and 2uj, with w/27r w 3 GHz, in a region in which the 
mw magnetic fields H[lo) and H{2uj) are maximal and 
parallel to each other. The w-niode of the cavity is fed 
by a train of mw pulses, with pulse width 5 /is, pulse 
repetition rate 200 Hz, maximum peak power ~ 50 W. 
The SH signal radiated by the sample is detected by a 
superheterodyne receiver. The cavity is placed between 
the poles of an electromagnet that generates DC magnetic 
fields up to w 10 kOe; two additional coils, externally 
fed, allow reducing the residual field within 0.1 Oe and 
working at low magnetic fields. All the measurements have 
been performed with Ho\\H{oj)\\H{2uj). The experimental 
apparatus is described in more details in Ref. 
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Fig. 1. SH signal intensity as a function of the temperature 
in the MgB2 sample, for //o = 10, 5, 2 Oe as displayed from 
the top. Input peak power ~ 30 dBm; the dashed line indi- 
cates the noise level that corresponds to « —75 dBm. The 
curves are normalized to the maximal intensity of the SH sig- 
nal detected at Ht^ = 10 Oe. The continuous lines have been 
obtained with the model discussed in the next section, using 
Ao/(So = 0.11, ao = 2 X 10-^Oe-\ oii/i = 8 x 10"^ 

Fig. 2 shows the SH signal intensity as a function of the 
DC magnetic field, at the two temperatures T = 37.6 K 
and T = 36.7 K. The results have been obtained in the 
ZFC sample on increasing for the first time. As ex- 
pected from symmetry considerations, the SH signal is 
zero at i7o = 0; on increasing the field it rapidly increases, 
exhibiting a wide maximum at fields of few Oe. Both the 
slope of the SH-vs.-iJo curves at Hq = and the position 
of the maximum depend on the temperature. 

By sweeping the magnetic field from zero up to a cer- 
tain value, Hmax, and back, the SH signal shows a hys- 
teretic behavior when Hmax reaches a threshold value that 
depends on temperature. Fig. 3 shows the field dependence 
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Fig. 2. SH signal intensity as a function of the static magnetic 
field, for T = 37.6 K and T — 36.7 K, as displayed from the 
top. Input peak power « 30 dBm. The curves are normalized to 
the maximal intensity of the SH signal detected at T = 37.6 K. 
The dashed line indicates the noise level. The continuous lines 
have been obtained by the model discussed in the next section, 
using the same values of the parameters of Fig. 1. 



In Fig. 4 we report the power dependence of the SH 
signal at Hq = 2 Oe, for two different values of the tem- 
perature. As one can see, the power dependence of the SH 
signal cannot be described by a mere n-order power law, 
indeed the slope of the SH-vs.-Pin line is not constant in 
the range of power levels investigated; in particular, at the 
investigated power levels the input-power dependence of 
the SH signal is less than quadratic. 

Input Peak Power (dBm) 




of the signal obtained sweeping Hq in the range ±15 Oe, 
for T = 37.6 K . The presence of a magnetic hysteresis, 
as well as the observation of a noticeable SH signal at 
Hq = 0, in the decreasing-field branch of the curve, sug- 
gest that trapped flux is present in the sample. So, we 
deduce that, even at temperatures near Tc, a critical state 
of the fluxon lattice develops. It is worth to remark that, 
on further increasing the DC magnetic field, the intensity 
of the SH signal decreases and, eventually, goes to zero 
when Hq reaches the value of the upper critical field. 
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Fig. 4. SH signal intensity as a function of the input peak 
power for Ho = 2 Oe and two different values of the temper- 
ature. The continuous lines have been obtained by the model 
discussed in the next section, using Ao/5o = 0.11, Qi=ao = 
2 X 10~^Oe~^. The dashed line indicates the noise level. 
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Fig. 3. SH signal intensity as a function of the static magnetic 
field, obtained sweeping Hq in the range ±15 Oe. Input peak 
power « 30 dBm. 



4 The Model 

The nonlinear electromagnetic response has been inves- 
tigated in conventional as well as in high- Tc supercon- 
ductors Several mechanisms responsible for the har- 
monic emission have been recognized; their effectiveness 
depends on the temperature, DC magnetic field and type 
of superconductor. At DC fields smaller than the lower 
critical field, the harmonic emission at T ^ Tc has been 
ascribed to nonlinear processes occurring in weak links [S] 
|ni[7] , while at temperatures close to Tc it has been ascribed 
to modulation of the order parameter induced by the mw 
magnetic field 'T','3','5',"^f3n' . 

Recently, we have investigated bulk MgB2 samples pre- 
pared by different methods and in all of them we have 
detected a near-T^ peak in the SH emission; however, in 
most of them the peak mingles with the low- T signal ■ 
As it has already been mentioned, the sample here investi- 
gated does not exhibit detectable SH emission at low fields 
and low temperatures, showing that the grain boundaries 
are not enough effective to give rise SH emission related to 
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processes occurring in weak links. This finding ensures the 
high sample quality and allows investigating the mecha- 
nism responsible for the harmonic emission near Tc in the 
MgB2 superconductor. 

Nonlinear emission at temperatures close to Tc is ex- 
pected in the framework of several models, some work in 
the weak-coupling limit pii:illl82j . others in the strong- 
coupling limit ,30^ . others in the framework of the two- 
fluid model ISElEl ■ In all of them, the harmonic emission 
originates from the time variation of the order parameter 
induced by the em field. 

From the models reported in Refs. pil3UU31U32) . the 
input power dependence of the nth-harmonic signal is ex- 
pected to follow an n-order power law. As it is shown 
in Fig. 3, we have observed a less than quadratic power 
dependence of the SH signal intensity; furthermore, the 
results reported in Ref. |25j show a less than cubic power 
dependence of the near- Tc third- harmonic signal in MgB2 . 
So, these models cannot justify the near- Tc harmonic emis- 
sion in MgB2. 

Results of second and third harmonic emission near 
Tc by high- Tc superconductors have been justified quite 
well in the framework of the two-fluid model, with the 
additional hypothesis that both the mw and the DC mag- 
netic fields, penetrating the sample in the surface layers, 
linearly perturb the partial concentrations of the normal 
and condensed fiuids 3 . Since the processes occur in the 
surface layer of the order of the field penetration depth, 
they are effective only at temperatures close to Tc, where 
the energy gap is not large and the penetration depth is of 
maximal extent. The model has been discussed for the first 
time by Ciccarello et al. 3 to account for the experimental 
data of third-harmonic emission by YBa2Cu307 (YBCO) 
single crystals in the absence of DC magnetic fields. Suc- 
cessively, it has been generalized to take into account the 
effect of the DC magnetic field and, further, justify the SH 
generation jJJBj . The model is valid in the limit of strong 
coupling, in which the approach of the two-fluid model 
can be reasonably used 13^1. It has been shown that the 
model accounts very well for the results obtained in high- 
Tc superconductors in the Meissner state and qualitatively 
those obtained in the mixed state 0. So, although it has 
been shown that the mechanism contributes to the har- 
monic emission also in the mixed state, others processes, 
related to the presence and motion of fluxons, may con- 
tribute to the nonlinear mw response of superconductors 
in the mixed state [MIT^ . 

Since the near- Tc SH signal radiated by MgB2 shows 
some peculiarities similar to those observed in YBCO and 
Bao.6Ko.4Bi03 (BKBO) crystals such as the input 

power dependence and the broadening of the peak on in- 
creasing Hq , one can reasonably hypothesize that in these 
superconductors the mechanism responsible for the near- 
Tc SH emission is the same. However, few differences have 
to be mentioned: the SH emission is more enhanced in 
YBCO than in MgB2 (at least one order of magnitude); 
the near- Tc peak is narrower in YBCO (~ Tc/30) than 
in MgB2 Tc/ 10); while in YBCO and BKBO the peak 
width, as well as the position of the maximum, depend 



on the input power, these effects have not been detected 
in MgB2. In the following of this section, we describe the 
model. 

Coffey and Clem |36j have elaborated a comprehensive 
theory for the electromagnetic response of type-II super- 
conductors in the framework of the two-fluid model of the 
superconductivity. It has been shown that the em field is 
characterized by a complex penetration depth 
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In Eqs. (2) and (3), Ao is the London penetration depth 
at T = 0; (5o = {c^ /2TTUjaoy/^ is the normal skin depth 
at T — Tc] Wo and 1 — wo arc the fractions of normal 
and superconducting electrons at _ffo = 0, in the Gorter 
and Casimir two-fluid model wq = (T/Tc)^. In Eq. (1), S-u 
is the effective complex skin depth arising from the vor- 
tex motion. At mw frequencies and at temperatures close 
to Tc it is reasonable to assume that vortices move in 
the flux-flow regime; in this case, using the expression of 
the viscous-drag coefficient given by Bardeen and Stephen 
123, it results Si = 25qBq/ Bc2(T), where Bq is the mag- 
netic induction due to the presence of fluxons. 

For a sample of thickness D much greater than the 
characteristic penetration depths, the mw magnetic in- 
duction averaged over the volume of the crystal is 



< B{t) > = 



D 



;3(A)costjt-f 3fi(A)sincji], (4) 



where Hi is the amplitude of the mw field. 

The Fourier coefficients of < B(t) > are given by 



1 



On = - 
TT 



TT 



2n 



< B{t) > cos{nujt)d{ujt), 



2tt 



< B{t) > sm{nujt)d{ujt). 



(5) 



(6) 



From a„ and 5„ one can calculate the induced magnetiza- 
tion oscillating at nuj and, eventually, the intensity of the 
nth-harmonic signal. 

The magnetic induction obtained using expressions (1)- 
(3) does not contain harmonic Fourier components be- 
cause all the equations involved up to now are linear. Non- 
linearity comes out on assuming that the mw field, pen- 
etrating the surface layers of the sample, modulates the 
partial concentrations of the normal and superconducting 
fluids 0. 

In the Meissner state, also the DC magnetic field pen- 
etrates in the surface layers of the sample and it can per- 
turb the partial concentration of the two fluids. On the 
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other hand, though we have performed measurements at 
low magnetic fields {Hq ^ 10 Oe), at temperatures very 
close to Tc the sample goes in the mixed state; in this 
case, a part of the DC magnetic field penetrates as flux- 
ons and the other is screened by the sample and is present 
in the surface layers. We suppose that both the mw and 
the screened DC magnetic field perturb the electron con- 
centrations in a similar way. Following this idea, it has 
been set 

w{t,Ho) = wo[l + \ao4TrM + aiHi coswt|], (7) 

where ao and ai are phenomenological parameters. 
By replacing wq in Eqs. ^ and © with w{t, Hq) of Eq. ((Tjl, 
the em induction inside the sample is expected to have 
harmonic Fourier components, which can be calculated by 
using Eqs. (O and |(HJ. The power emitted by the sample 
at nuj is expected to be: 

Pnc. oc al + bl. (8) 



5 Discussion 

By using Eqs. H4I8(I we have deduced the expected SH 
signal intensity. The integrals of Eqs. (0) and © cannot 
be solved analytically; they have been solved numerically 
with the proviso that w{t,HQ) — 1 whenever, because of 
the contribution of the perturbative terms in Eq. Q, its 
instantaneous value becomes greater than one. This con- 
dition is necessary since the concentration of normal elec- 
trons cannot exceed one and no modulation must occur 
when all electrons are in the normal state. 

The calculations have been carried out in the follow- 
ing approximation: |47rM| equal to Hq for Hq < Hci{T); 
\A7tM\ linearly decreasing from iJci (2^) to zero for iJci (T) < 
Ho < Hc2{T). 

The expected results depend on the ratio Ao/<^o, the 
temperature dependence of wq and the parameters ao and 
ai. In particular, Xo/So and the temperature dependence 
of Wo determine the peak width, ao and ai the field and 
power dependence of the SH signal, respectively. Since the 
width of the near-Tc peak is most likely affected by the 
inhomogeneous broadening of the superconducting transi- 
tion, we have calculated the averaged em induction (Eq.^ 
taking into account the distribution of Tc over the sample. 
By using reasonable values of Aq /So and the temperature 
dependence of wo expected from the Gorter and Casimir 
two-fluid model, we have obtained a near-Tc peak much 
narrower than the one experimentally observed. On the 
other hand, different authors |19ll20ll38ll39j have shown 
that the temperature dependence of the field penetra- 
tion depth in MgB2 cannot be accounted for by cither the 
Gorter and Casimir two-fluid model or the standard BCS 
theory. A linear temperature dependence of the condensed 
fluid density, in a wide range of temperature below Tc , has 
been reported, which has been justified in the framework 
of two-gap models for the MgB2 superconductor [T??li^ . 

Prompted by these considerations, we have calculated 
the averaged em induction over the sample, from Eq. 



supposing a linear temperature dependence of wq and 
using Xq/So, ao and ai as parameters. The continuous 
lines of Figs. 1 and 2 show the temperature and field de- 
pendence of the SH signal, respectively, expected from 
the model. The best-fit curves have been obtained with 
Xo/5o = 0.11, ao = 2 X lO^^ Oe"! and aiHi = 8 x lO-^. 
Since the amplitude of the mw magnetic field at which we 
have performed the measurements is about 4 Oe, we may 
estimate ai to be of the same order of ao. By using the 
same value of Xo/Sq and setting ai = ao = 2 x 10"'^ Oe~^ 
we have obtained the continuous lines of Fig. 4 for the 
power dependence of the SH signal intensity. As one can 
see, the results expected from the model, hypothesizing 
that the normal and condensed fluid densities linearly de- 
pend on temperature, satisfactorily agree with the exper- 
imental data. We would remark that, by using the same 
values of the parameters before mentioned, we have found 
that, contrary to what occurs in other superconductors 
EEEI, the width of the SH near-Tc peak in MgBs is not 
appreciably affected by the amplitude of the microwave 
magnetic field. This finding, which agrees with the ex- 
perimental observations, may be ascribed to the different 
values of the parameters used to fit the experimental data 
of MgB2, with respect to those obtained in other super- 
conductors. 

Since Ao /So — ycijr/2, by supposing Ao of the order of 
100 nm wc may estimate the scattering time of the normal 
electron to be r ~ 10~^^ s, which is of the same order of 
that reported in the literature for MgB2 001 ■ However, we 
would remark that we have assumed a linear temperature 
dependence of the condensed fluid density, which could 
be no longer valid at low temperatures; so, it is possible 
that the value of Xo/Sq we found is not exactly equal to 

The hysteretic behavior shown in Fig. O suggests that 
mechanisms related to trapped magnetic flux are effective. 
On the other hand, it is reasonable to hypothesize that at 
temperatures near Tc the sample goes in the mixed state 
for magnetic flelds of the order of 10 Oe. The presence of 
the magnetic hysteresis suggests that a critical state of the 
fluxon lattice develops in the sample up to temperatures 
near Tc. Since the model here reported does not consider 
such effects, we did not fit these results. When the sam- 
ple is in the critical state, the magnetic induction, Bq, is 
different for applied magnetic field reached at increasing 
and decreasing field, causing a magnetic hysteresis in the 
SH signal. The sharp minima observed after the inversion 
of the field-sweep direction may be ascribed to the change 
of the critical current in the sample skin layer |30|. 



6 Conclusions 

In conclusion, we have reported a set of experimental re- 
sults on second harmonic generation at temperatures close 
to Tc in a bulk sample of MgB2, prepared by the reac- 
tive liquid Mg infiltration technology. The SH signal has 
been investigated as a function of temperature and ampli- 
tudes of the DC and mw magnetic fields. The results have 
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been discussed in the framework of a phenomenological 
model previously elaborated to account for the nonlin- 
ear emission near Tc by YBCO and BKBO crystals. The 
model assumes that at temperatures close to Tc the mw 
and the DC fields, penetrating in the surface layers of the 
sample, perturb the partial concentrations of the normal 
and superconducting electrons. We have shown that, in 
order to account for the experimental results, it is essen- 
tial to suppose that in MgB2 the densities of the normal 
and condensed fluids linearly depend on the temperature. 
Such finding agrees with experimental results of the tem- 
perature dependence of the penetration depth reported 
in the literature, which have been justified in the frame- 
work of two-gap models for MgB2 superconductor. We 
have shown that the enhanced nonlinear effects detected 
near Tc in MgB2 have the same origin of those observed in 
YBCO and BKBO crystals. Though the superconductors 
have different nature, the description of the electromag- 
netic properties in terms of the two-fluid model seams to 
be appropriate. The different peculiarities of the near- Tc 
peak in the SH response, among the investigated super- 
conductors, have to be ascribed to a different temperature 
dependence of the normal and condensed fluid densities. 
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